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Abstract.  The  existence  of  weak  solutions  to  the  equations  proposed  by  Charney 
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1.  Introduction.  In  this  paper,  we  examine  mathematical  properties  of  an  equation 
arising  in  the  theory  of  ocean  circulation.  In  order  to  understand  the  role  of  this  problem 
in  oceanography,  a  brief  review  of  the  subject  is  necessary.  The  first  successful  attempt  to 
provide  a  mathematical  model  of  the  mid-latitude  ocean  currents  was  made  by  Stommel 
[11]  in  1948.  He  showed  conclusively  that  a  Gulf-Stream-like  intensification  on  the  western 
side  of  an  ocean  basin  could  be  explained  by  the  so-called  0— effect.  This  is  the  geophysical 
terminology  for  the  latitudinal  variation  of  the  normal  component  of  the  Earth’s  rotation. 
Aside  from  this  variable  Coriolis  force,  the  other  forces  which  entered  into  Stommel’s  model 
were  those  due  to  the  pressure  gradient,  the  surface  winds  and  friction.  This  last  force 
was  taken  to  be  proportional  to  the  velocity  fields.  All  the  effects  of  density  stratification 
were  neglected  by  making  the  assumption  that  the  ocean  was  homogeneous.  Finally,  by 
working  with  vertical  averages,  Stommel  essentially  treated  the  ocean  circulation  as  a  two- 
dimensional  horizontal  motion.  Somewhat  surprisingly,  Stommel’s  ad  hoc  linear  model  was 
shown  later  to  provide  an  accurate  description  of  an  experimental  set-up  [10]. 


The  subsequent  work  in  the  field  has  attempted  to  overcome  the  two  oversimplifica¬ 
tions  of  Stommel’s  model.  Namely,  to  take  into  account  inertial  forces  which  are  known 
to  be  important  in  Gulf-Stream  like  boundary  layers  and  to  represent  more  adequately 
the  complex  mixing/dissipative  processes.  These  two  broad  generalizations  were  initiated 
by  Charney  [4]  and  Munk  &  Carrier  [9]  respectively.  Although  these  early  papers  make 
use  of  analytical  techniques,  the  bulk  of  the  recent  work  in  this  field  has  relied  on  numeri¬ 
cal  techiques  to  study  the  nonlinear  problems  arising  in  the  mathematical  formulation  of 
wind-driven  ocean  circulation  models.  The  papers  by  Bryan  [3]  and  Veronis  [13], [14]  fall 
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within  this  category. 


In  this  paper,  we  study  the  Charney’s  extension  of  Stommel’s  model,  namely 


(1.1) 


cA  ip  + 


dip 

dxi 


+  *( 


dip  dAip 
dxi  dx 2 


dip  d  a  ip 

dx 2  dx 1 


in  a  bounded  domain  Cl  in  R2  with  Dirichlet  boundary  conditions.  We  prove  existence 
of  weak  solutions  to  (1.1).  These  weak  solutions  are  obtained  as  limits  of  solutions  of 
auxiliary  equations  with  artificial  viscosity  and  artificial  boundary  conditions.  Physically, 
this  auxiliary  problem  corresponds  to  the  addition  of  side  wall  friction  and  stress-free 
boundary  conditions.  We  derive  uniform  L°°  bounds  for  the  solutions  of  the  auxiliary 
equations.  This  procedure  cannot  yield  classical  solutions  because  the  artificial  boundary 
conditions  give  rise  to  boundary  layers  in  the  classical  limit.  Nevertheless,  we  expect  the 
solutions  to  be  classical  in  certain  parameter  ranges. 


The  same  method  was  used  by  Yudovitch  [15]  in  his  paper  dealing  with  the  time 
dependent  two-dimensional  Euler  equations.  This  method  cannot  be  extended  to  the  time 
independent  Euler  equations.  Actually,  it  is  known  [8]  that  the  time  independent  Euler 
equations  have  no  solutions  for  two  and  three  dimensional  axisymmetric  domains  and 
generic  driving  forces.  However,  because  of  the  presence  of  the  bottom  friction  term  eA  ip, 
the  method  can  be  used  for  equation  (1.1). 


Once  the  estimates  for  the  stationary  problem  (1.1)  are  understood,  the  time  depen¬ 
dent  version  of  this  problem 


(1.2) 


dAip 

di 


A  /  ,  w  ,  D  dip  dAlp 


dzi  dx 2  dx 2  dx\ 

with  given  initial  conditions  on  A  ip  can  be  obtained  by  applying  in  a  straightforwardor 
manner  the  method  of  Yudovitch.  The  /3— effect  term  dip /dx i  does  not  create  any  serious 
difficulties.  One  can  show  that  (1.2)  has  unique  global  solutions. 


The  problem  of  describing  the  set  of  stationary  solutions  is  still  open.  We  expect  the 

solution  to  be  unique  only  when  e  is  large  compared  to  R  \  ,  '  '  ~  - - H 

or  .Distribution/ 


□  □ 


2.  Preliminaries.  Let  ft  be  an  open  bounded  set  of  R2  with  sufficiently  smooth 
boundary  0ft.  We  consider  the  system 


dip  du  dip  du>  .  '  O 

dx\  dl2  dl2  dx\ 

Aip  —  oj  in  ft 

ip  =  0  on  0ft 

where  f(x),c  >  0  and  R  >  0  are  given. 

In  this  paper  we  study  the  existence  of  weak  solutions  of  (2.1). 

We  denote  by  if* (ft)  the  usual  Sobolev  spaces  of  order  s  and  by  Hq(Q)  the  closure  of 
<7o°(ft)  in  the  H1((i)  norm.  We  define 


A  =  -A 

to  be  the  negative  Laplacian  with  domain  D(A)  =  iT2(ft)  H  J3o(ft).  It  is  well  known 
that  A-1  is  a  compact  linear  self-adjoint  positive  operator  in  L2(ft)  (cf.  [2],  [7]).  The 

spectrum  of  A  consists  of  an  infinite  sequence  0  <  Ai  <  A2  <  . . .  of  eigenvalues  counted 
according  to  their  multiplicities;  An  — »  oo  as  n  — ♦  oo;  the  eigenfunctions  {u>»}  provide  an 
orthonormal  basis  in  Z<2(ft).  Finally,  there  exists  a  scale  invariant  constant  cq  snch  that 
|ft|  =  CflAj-1  where  |ft|  denotes  the  area  of  ft.  The  scalar  product  and  norm  in  L2(Sl)  are 
denoted  by  (•,  •)  and  |  •  |  respectively.  The  scalar  product  in  /^(ft),  in  view  of  Poincare’s 
inequality,  is 


((u,v))=  f  Vt i(x)  Vv(x)dx  Vu,u€Fo(ft), 

Ja 

and  the  corresponding  norm  is  denoted  by  ||  •  ||.  The  norm  in  Lp{il)  for  1  <  p  <  oo  is 
denoted  by  ||  •  Up.  It  is  known  (cf  [2])  that  on  D(A)  the  if  2(ft)  norm  is  equivalent  to  the 
| A  •  |  norm,  i.e.  there  exists  a  constant  cj  >  0  snch  that 


<  ||u||h*(o)  <  ci|i4u|  Vu  e  D{A). 


Moreover,  D(A1!2)  =  Hq( Q)  and  ||  •  ||  =  l^1^2  •  |. 
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The  bilinear  bounded  operator  J  :  H *(0)  x  H1( fl)  — »  ^(fl)  is  defined  as 


dip  du  dip  d(jJ  i 

=  T- -ir-  ir  v  r,"  £#'<«)• 

C/2l  1712  t7X2  vf*l 


We  need  the  following  inequalities  which  correspond  to  various  continuity  properties 
of  the  operator  J. 


Proposition  2.1.  Let  si,  32,  *3  >  0  satisfy 


si  +  *2  +  *3  >1  if  a,  1 


for  all  i  =  1, 2, 3 


31  +  32  +  ^3  >1  if  «»  =  1  for  some  i  =  1,2,3. 


Then,  there  exists  a  scale  invariant  constant  c(si,  32**3)  such  that 


(2.3)  KW,  «),»)|  < sj)x(,“+'I+“),2!!^|!ff>1  +  i  ■  |H!a„  +  j  •  IMU.S(0). 


The  reader  is  refered  to  [12]  and  [5]  for  the  idea  of  the  proof.  We  also  note  that 


w),v)  =  ~(J(u,ip),v) 


for  every  1 1>  G  ff<1  +  1(fi).w  G  +  and  v  G 


Proposition  2.2.  For  every  ip  G  H2(Q),  u  G  JEf  1(S1)  and  v  G  Hq(U)  we  have 


( J(lp,Uj),V )  =  —(J(ll>,v),  U>). 


Proof:  Because  of  proposition  2.1  '  a  sufficient  to  show  that  (2.5)  holds  for  v  €  Co°(fi) 
and  ij>,u  €  Cao(n).  Let  tT  be  the  vectc.  [~dij)/dx2,  dr/>/dx\);  then  one  can  write 


u/),v)  =  f  V  ■  (u(x)u{x))  v(x)  dx 

Ja 

=  —  I  (5(i)-Vv(i))a >(x)dx  . 

Ja 

The  following  corollary  is  an  immediate  consequence  of  the  above  proposition. 


Corollary  2.3. 

(2.6)  (J(Vs w),w)  =  0  Vtf  €  tf2(fi),u>  €  H^{Q) 

(2.7)  (J(V>, w)Jc^)  =  0  W  e  H2(Q),u  e  D(A),  p  =  1,2, . 


Proof:  (2.6)  is  a  direct  consequence  of  (2.5).  Since  iT2(Q)  is  a  Banach  algebra  (cf. 
[1])  one  can  easily  verify  that  Jp  €  D(A),  for  p  =  1,2,...,  whenever  u  €  D{A).  We  can 
then  use  (2.5)  to  deduce  that 


(/(Vsu),!!^)  =  -(JT(^,u^),u;), 
while  from  direct  computations 


w)  =  -p(J(V>,u>),u;),  p=  1,2,... 

which  verifies  (2.7),  in  view  of  (2.6).  Let  us  note  that  (2.6)  also  holds  for  i{>  €  C1(Q), 
u>  £  C1  ( 0)  and  normal  to  dCl. 


3.  A  stationary  problem  with  artificial  viscosity.  In  this  section  we  consider  a 
singular  perturbation  of  the  stationary  problem  (2.1)  obtained  by  adding  an  artififial  vis¬ 
cosity  term  -vAw  and  imposing  homogeneous  Dirichlet  boundary  conditions  on  u.  The 
artificial  viscosity  equation  provides  an  approximate  solution.  Uniform  bounds,  indepen¬ 
dent  of  the  viscosity  v  will  enable  us  to  pass  to  the  weak  limit.  In  this  process  the  boundary 
condition  on  u;  is  lost.  The  auxiliary  problem  is: 


(3.1a) 

-v Aa>  +  +  dil>/dx\  -f  RJ(%I),u)  =  / 

in  0 

(3.16) 

—  +  u>  =  0 

in  Q 

(3.1c) 

■0  =  0,  U>  =  0 

on  #Q, 

where  v 

is  given  such  that  0  <  v  <  e3/8  and  /  G 

3.1  Existence  of  solutions  of  the  perturbed  problem.  We  restate  the  problem 
thus:  find  (u;,V>)  €  D(A)  x  D(A)  such  that 


(3.2a)  vAu  +  eu}  +  dij)/dx\  +  RJ(t/’,w)  =  f, 

(3.2  6)  Aij>  +  u>  =  0. 

It  is  obvious  that  every  solution  (w,  rj’)  of  (3.2)  is  a  solution  to  (3.1)  in  the  distribution 
sense. 

The  following  a  priori  estimates  of  the  solutions  of  (3.2)  will  be  needed. 


Lemma  3.1.  Let  (w,ip)  €  D{  A)  x  D(A )  be  a  solution  of  (3.2).  Then 


2 


e 


and 


(3.6)  v\Au[  <  2|/|  +  2V2!/iA-11',(<,A1)[l  +  °)C’  l/l’tf.M.)], 

where 

«Tl(«,A,)=  (1  +  -2_). 


Proof:  Let  (w,^)  £  D{A)  x  D(A )  be  a  solution  to  (3.2).  Taking  the  scalar  product  of 
(3.2a)  with  if)  we  get: 


(3.7) 


v{Au>,  i>)  +  e(o;,  i>)  +  (~,  i>)  +  R{J{tP,  ui),  i>)  =  (/,  V>)- 


In  view  of  (2.4)  and  (2.5)  and  of  the  fact  that  (difr/dx i,  i/>)  =  0,  (3.7)  reduces  to 


or  because  of  (3.2b) 


v(Au,  ij>)  +  e{ u,  if>)  =  (/,  V'), 


e||^||2  <  ,|uf  +  I/!  1'01- 


Using  first  Poincare’s  inequality  we  can  write 


e\\i>\\2  <  v\uj\2  +  )/)  X1  1/2  jlV’H 

then  by  Young’s  inequity 


and  finally  since  0  <  u  <  c3/S,  we  deduce 


||^||2<(2I//e)H2  + 


l/l2 


e2Ai 


If  we  take  the  scalar  product  of  (3.2a)  with  w  and  use  (2.6)  we  obtain 


u\\u\\2  +  c|w|2  =  -{dip/dxi,  w)  +  (/,  w) 


and  by  Young’s  inequality 


<  11-011  H  +  l/l  M, 


W  +  \J? 


From  the  above  inequality  and  (3.8)  we  deduce  (3.3).  From  (3.3)  and  (3.8)  we  can 
obtain  easily  (3.4).  Finally,  (3.2b)  and  (3.3)  directly  imply  (3.5). 


To  establish  (3.6)  we  take  the  scalar  product  of  (3.2a)  with  Au >,  from  which  we  see 

that 


v\Au\2  +  e|H|2  <  \(dif>/dzi,  Aw)]  +  Aw)|  +  | (/,  Aw)|, 

<  ill’ll  |Aw|  4-  ii|(J(^,w),Aw)|  +  |(/,  Aw)|. 


By  using  (2.3)  for  si  =  32  =  1/2  and  33  =  0,  we  get 


t/|Aw]2  +  e||w|l2  <  (|!V»!i  +  Rc(  1/2, 1/2,0)  IM|jj*/*(n)  +  1/1)1^ 


and  by  an  interpolation  inequality 


HAuf  +«MI  <  {H\\  +  *c(  1/2,  l/2,0)cM|''» +  l/l)|A*l. 

and  hence 


^|Au>|2  +e|ju;jj2  <  (|M|  + 


^U/2.  !/2. 0)c2  iiV'll  \M\\M\  +  I/I)  \Au\. 


Making  use  of  (3.3), (3.4)  and  (3.5)  yields  (3.6). 


Lemma  3.2.  There  exists  at  least  one  solution  to  problem  (3.2).  Moreover,  every  solution 
of  (3.2)  satisfies  (3.3)-(3.6). 


Proof:  One  approach  is  to  use  the  Galerkin  approximation  method  based  on  the 
eigenfunctions  of  the  operator  A  as  in  Constantin  fc  Foias  [5]  and  Temam  [12]  to  show  the 
existence  of  a  solution  to  (3-2).  The  crucial  point  in  this  approach  is  to  establish  similar  a 
priori  estimates  to  the  ones  in  (3.5)-(3.6)  for  the  approximate  solution. 


Another  approach  is  by  the  Leray-Schauder  degree  theory.  Indeed,  problem  (3.2)  is 
equivalent  to 


<V[V 

i/~l[eu>  +  dtp/dxx  +  R 


where 


(v  [eu -t  oy/oxi -t  njy\p,w)]\ 
R(y>,u>)={  j. 

By  Rellich’s  Lemma  and  (2.3),  one  can  show  that  the  nonlinear  mapping 

K  :  D(A)  x  D(A)  -  L2{Q)  x  L2(0) 


is  compact.  Therefore,  making  use  of  the  a  priori  estimates  (3.5)-(3.6),  we  can  conclude 
that  (3.9)  has  at  least  one  solution. 

Remark  3.3:  If  /  €  then  every  solution  G  D(A)  x  £>(.4)  of  (3.2) 

satisfies  (3.1)  in  the  classical  sense,  namely  G  C°°(fi). 

3.2  Uniform  L°°  bounds  for  the  artificial  viscosity  problem  In  (3.3),  (3.4)  and 
(3.5)  we  gave  estimates  for  |u;|,  \\if>\\  and  j  Atp\  which  were  independent  of  v.  In  this  section 
we  shall  derive  uniform  (i.e.  independent  of  v)  L°°  estimates  for  every  solution  of  (3.2) 
and  for  0  <  u  <  e3/8.  We  recall  first  the  following  known  result  in  potential  theory  (see 
e.g.  [6] ). 


Theorem  3.4.  Let  G(x,  y)  denote  the  Green  function  of  the  Laplacian  operator  in  the  do¬ 
main  with  Dirichlet  boundary  conditions.  Then,  there  exists  two  scale  invariant  constants 
c 3,  c4  such  that: 


(3.10) 
and 

(3.11) 


\G(xfy)\  <  c3  (1  +  |log(Aj/2|z  -  y\)\)  Vz,y  G  Q,x  ^  y , 


Vfc  =  1,2,  x,yen,x?y. 


dxi 


ix  -  v ! 


Moreover  since 


it  follows  that 
(3.12) 


4>(x)  =  f 
d<t>(x)  _  f 

Jn 


G{x,y)Ac}>(y)dy , 


dG(x,y) 

dxh 


A<p(y)dy. 


Theorem  3.5  (Uniform  £°°  bounds).  Let  (u>,V’)  G  D{A)  x  D(A)  be  a  solution  of  (3.2) 
(or  equivalently  a  solution  of  (3.1).  Then: 


an 


(3.14) 


HVV'lloo  <  4 v^rc4cj/4  \x  1/4  v/iHUH. 


Note  that  (3.3)  provides  a  uniform  upper  bound  for  |u>| . 


Proof:  Let  6  >  0.  We  denote  by  B((z)  the  ball  in  RJ  that  is  centered  at  z  with  radius 
6.  Prom  (3.1a)  and  (3.13)  we  have 


Vi>{x)=  f  VxG(x,y)u>(y)dy, 
Jn 


which  on  account  of  (3.11)  implies 


|v*(*)|<2c  / 

J n  I®  -  y\ 

r  f  M»)l  . .  ,  f 


! VV’(*)|  <  2c4 


I  f  Ms 

*■ J(ir\B,(z )  I*  - 


ldy+f  M?>1 

y\  Jn\B,(x)\x  -  y\ 


^4 


Using  Holder’s  inequality,  we  see  that  for  odd  integers  p  >  3 


|V^(*)|<2c4{|H|F+1(jr  \*-y\  {p+l),pdy)r/p+x  +  ^|fl|1/2}, 


|V^(*)|  <  2c4{||a;||J>+I(2x)^+1^-l^+1)  +  ^|fi|1/2}. 


and  therefore 


(3.15) 


IIVVIU  <  (4»c.)|H|,+1  s>-'i’+l  +  2c4!m1'Jy. 


In  order  to  get  an  estimate  for  ||u;||j,+i,  we  multiply  (3.1a)  by  u^(x)  and  integrate  over  0. 
This  leads  to 


—  uf  Au:(x)iJ!,(x)dx  +  e  f  u?+l(x)dx  -f  R  (J(ifi,u;),u>p)  =  —  (dip/ 8x1,0?)  +  (/,w) 

Jn  Jn 


or,  after  integrating  by  parts  and  using  (2.7),  to 


f  |  Vu;(z)|2  uP  l(x)dx- f-c/  <jjp+x(x)dx  +  R  =  -(dy>/dxu  u>p)  +  (f,u>). 


i, rp*  *.«  * '  ■  < 


Using  Holder’s  inequality  as  well  as  the  (act  that  p  is  an  odd  integer,  we  conclude  that 


‘IKS  < 


\\  w  T  ||< 


or 

(3.16)  e|!«!i;:i  <  (l|v*||„|n|‘"+1  +  II/IU0IH5+.- 

Substituting  (3.15)  in  (3.16)  we  see  that 


(3.17) 


€  UM 


ip+1 


<  ||/||,+1  -I-  (41rc4)|n|1/(?4l)«5(y-1)/(3’+1)||u;||J)+1+  2c4 


If  we  choose  6  =  (e/8irc4)|n|  l^p  then  (3.17)  implies  that 

2  4XC4  ^8lTC4  jQjl/2|P+1/»,_1|Q|l/p+lj(^|^ 


(3.19) 


Mlp  +  i  < 


p+i 


e  •  e  e 

Passing  to  the  limit  as  p  — *  oo,  we  conclude  that 


(3.19) 


2  ..  Z2rc\  .  ,/a 

IHl»  <  -  Il/Iloo  +  M  "M- 


After  replacing  |0 j  in  the  above  expression  by  cq  1  and  using  (3.3),  we  obtain  (3.13) 


In  order  to  derive  (3.14),  we  observe  that  (3.15)  holds  for  every  6  >  0.  Therefore 

(3.20)  || VV^iioo  <  (4xc4)  Halloo  6  +  2c4  |ft|1/2  !jL 

Minimizing  the  right  hand  side,  we  deduce  that 

IIVVHU  <  4v/2tc4  in^4  v/iHUH, 

and  after  replacing  |fl|  by  c©  x,  we  obtain  (3.14). 


4.  Existence  of  weak  solutions.  Let  {i'fc}fcL1  be  an  arbitrary  sequence  of  real 
numbers,  i/*  €  (0,  c3/8)  for  k  =  1. 2, . . and  such  that  —*  0  as  k  — *  oo.  Let  (a>k,  i/>j,)  be  a 
sequence  of  solutions  to  problem  (3.1)  corresponding  to  v  =  Uk  for  k  =  1, 2, . . .  respectively. 
By  Lemma  3.2,  we  know  that  such  sequences  exist. 


Theorem  4.1.  There  exists  a  ip  €  D(A)  and  an  w  £  that  are  weal  solutions  of 

the  stationary  problem  (2.1),  i.e. 


(4.1) 


(  fu>  +  dip/dx  1  -f  /?  V  •  (wu)  =  /  in  0 

=  ^  in  ft 

■0  =  0  ob  5ft 


w here  u  =  (  —  dip  j  dxi.  dip/dx\).  Furthermore,  ip  and  u;  satisfy: 

2|/l  t'l/2/ 


(4.2) 

(4.3) 

(4.4) 


M  < 


AV^e.Aj), 


2  fi4  1/2 

ll^lloo  <  -  ll/llcc  +  ^  |/i  K\>\e,\x), 

!v^|U<4v^c4cJ/4A-1/4  v/jHU  H. 


Proof:  Let  (u)it,i>h)  be  as  mentioned  above.  From  (3.5)  we  know  that 

\Ai>k\  <  ~  Kl/2(t,  Ai). 

Therefore,  there  exists  a  subsequence,  say  {V’fcilt  which  converges  weakly  in  ff2(ft)  to 
ip  €  H2(Q),  and  by  virtue  of  Rellich’s  Lemma,  ipitl  converges  strongly  to  ip  £  and 

hence  ip  G  D(A).  From  (3.13)  we  know  that 


IK, II-  <  f  ll/IU  +  I/I 

and  therefore 

(4.5)  K,n,<mi,/'  (?„/,. +  vp = 


By  means  of  (4.5)  we  can  inductively  find,  for  every  p  =  2, 3, ...  a  subsequence  {u^} 
of  {wfc}  which  converges  weakly  in  ZA^ft)  to  the  same  limit  u>  €  IA’( ft)  for  every  p  =  2,  3, . . .. 
In  particular,  we  have 


|MI„  <  liminf  ||u;fc  ||  Vp  =  2,3,... 

«  p— »oo 


which  by  (4.5)  entails  that 


(4.6)  imi,  <  mi 


Vp=  2,3, 


Passing  to  the  limit,  we  deduce  (4.3).  The  derivations  of  (4.2)  and  (4.4)  follow  simply 
from  (3.3)  and  (3.14).  We  shall  omit  their  proofs. 


In  conclusion,  the  diagonal  subsequence  }  converges  weakly  to  u  €  £**(0)  for 
every  p  =  2,3,...,  and  {i>kh}  converges  weakly  to  G  H2((l)  and  strongly  in  Hl(Q). 
Because  converges  strongly  in  £2(Q)  to  V^>  and  converges  weakly  in  £2(fl),  it 

follows  that  the  product  (VV,fck)(w*lk)  converges  to  (VV>)(w)  in  distribution  sense.  Since 
•{Vk4,V,kt}  solves  (3.1)  for  v  =  i/*,  by  passing  to  the  limit  we  can  verify  that  {v,tp)  solves 
(4.1). 


Up  to  now,  the  parameter  R  did  not  enter  at  all  in  our  estimates.  In  the  next  theorem, 
we  give  an  upper  bound  to  the  “diameter”  of  the  set  of  stationary  solutions  of  problem 
4.1. 


Theorem  4.3.  Let  (V>  1 ,  wi )  and  (V>2,  a>2)  be  two  weak  solutions  to  problem  (4.1)  satisfying 
(4.2)-(4.4).  Then 

(4.7)  lh>,->MI<^l/|,K1(€,A1)  (gll/IL+^f  I/I  *;'*(<■  A,)j. 


Proof:  From  (4.1)  we  have: 


eA(V>i  -  ^2)  +  ^-(^1  ~  ^2)  +  i2div(u;1ui  -  u>2u 2)  =  0, 
or 

Q 

eA(-0i  -  ^2)  +  ^-(V’l  “  i>i)  +  i?div(w1(it1  -  u2)  +  (o>i  -  u;2)tf2)  =  0. 

We  form  the  scalar  product  with  (^1  -  ^2)  and,  since  V’i  —  V>2  G  i7g(Q),  we  can  integrate 
by  parts  and  obtain: 


ellV’i  -  +  #(£2(^1  -  U72 ),  V(V’i  -  i>2))  =  0. 


From  (4.8)  we  get: 

ellV’i  -  V’slt2  <  R  ||«i  -  wjIIoo  |u3|  ||-01  -  0a||, 

and  hence 

(4.9)  ||0i  -  0a II  <  ~  ||«i  -  wjHoo  ||0ail- 

As  a  consequence  of  Lemma  3.1  and  Theorem  4.1, 

Will  <2  I/I’  *!(<,*,) 

Substituting  this  expression  and  (4.3)  in  (4.9)  we  arrive  at  (4.7). 
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